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We investigate the entanglement dynamics of continuous- variable quantum channels in terms of 
an entangled squeezed state of two cavity fields in a general non-Markovian environment. Using 
the Feynman- Vernon influence functional theory in the coherent-state representation, we derive an 
exact master equation with time-dependent coefficients reflecting the non-Markovian influence of 
the environment. The influence of environments with different spectral densities, e.g., Ohmic, sub- 
Ohmic, and super-Ohmic, is numerically studied. The non-Markovian process shows its remarkable 
influences on the entanglement dynamics due to the sensitive time-dependence of the dissipation and 
noise functions within the typical time scale of the environment. The Ohmic environment shows a 
weak dissipation-noise effect on the entanglement dynamics, while the sub-Ohmic and super-Ohmic 
environments induce much more severe noise. In particular, the memory of the system interacting 
with the environment contributes a strong decoherence effect to the entanglement dynamics in the 
super-Ohmic case. 

PACS numbers: 03.65.Yz, 03.67.-a 



I. INTRODUCTION 

Quantum teleportation incorporating the classical 
communication theory with a unique characteristic of 
quantum mechanics, quantum entanglement, has re- 
ceived tremendous attention in the study of quantum 
communication in the past decade [J [2j In quan- 
tum teleportation protocols, a necessary ingredient is the 
quantum channel, which is realized through an entan- 
gled quantum state of two systems separated between 
the sender and the receiver. Theoretically, both the 
discrete, (for example, two polarized photons, two-level 
atoms, orthe spins of electrons etc.)[l, 0, 0, 0, H, Q 
and the continuous-variable (coherent and squeezed op- 
tical fields) [13, [HI entangled states are equally useful 
for a quantum channel. Practically, compared with the 
discrete- variable entangled state, the continuous- variable 
entangled state may be more efficient because it has 
less decoherence [l^, [13, [3 • Continuous- variable entan- 
gled states can be traced back to the original paper on 
quantum entanglement by Einstein, Podolsky, and Rosen 
[l5| . where entangled states of the common eigenstate 
of relative position and total momentum of two parti- 
cles were proposed. Such ideal entangled states can ac- 
tually be realized by a two-mode squeezed state of op- 
tical fields in the large limit of the squeezing parame- 
ter. In fact, the entangled two- mode optical squeezed 
state has been successfully produced via the nonlinear 
process of parametric down conversio n I l6l| . This trig- 
gered a variety of experiments [H, [H, Il4j applying such 
an entangled state to quantum teleportation. The en- 
tangled two-mode optical squeezed state has been of 
key importance as an entangled resource for practical 
implementations of quantum-information protocols 
However, a realistic analysis of any quantum channel 



must take into account the noise effect from its environ- 
ment. There has been an increasing interest in describing 
continuous-variable entan glem ent d y namics under noise 

[13, M [H, Hi, [m, [H, mm, ^m- 

The traditional approach to studying the environment- 
induced noise effects treats the interaction between the 
quantum system and its environment perturbatively, 
which yields approximate equations of motion such as 
Rcdfield or master equations under the Born-Markov ap- 
proximation [23, [H, [13|- Although this treatment has 
been widely employed in the field of quantum optics, 
where the characteristic time of the environmental corre- 
lation function is much shorter compared with that of the 
system investigated [l^l , its validity is experiencing more 
and more challenges in facing new experimental evidences 
[30| . Moreover, the Born-Markov approximation is in 
general invalid in dealing with most condensed-matter 
problems, for example, a quantum system hosted in a 
nanostructured environment 31, 32, 33, 34, 35], because 
large coupling constants and long correlation time scales 
of the environment both require a non-perturbative de- 
scription. Therefore a nonperturbative description of the 
non-Markovian dynamics in open quantum systems has 
attracted much attention over recent years [36|. 

In fact, recently, non-Markovian processes have been 
extensively studied in the entanglement dynamics of two 
continuous-variable systems, such as two harmonic os- 
cillators or two-mode electromagnetic fields, interacting 
with bosonic environments [H, [2I [H [H, [Si] . In [l^, 
the non-Markovian entanglement dynamic of two-mode 
Gaussian states is studied based on the master equa- 
tion derived perturbatively using the projection opera- 
tor method up to the second order with respect to the 
system-reservoir coupling constant (which actually cor- 
responds to the Born approximation). In (23[ the dy- 
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namics of two harmonic oscillators interacting with two 
uncorrelated reservoirs was formulated based on the Hu- 
Paz-Zhang master equation of quantum Brownian mo- 
tion [39| but the non-Markovian entanglement dynam- 
ics is analyzed up to the second order of the system- 
reservoir coupling constant with Ohmic reservoirs. More 
recently, an exact master equation for two coupled har- 
monic oscillators linearly interacting with a common 
reservoir has been derived usin g th e Feynman- Vernon 
influence functional theory [stI. l38l. [soj where the de- 
coherence and disentangled dynamics of a bipartite dis- 
placed Gaussian states is studied within the Markovian 
approximation 24 1 . In a very recent paper [26| , the non- 
Markovian master equation of two coupled harmonic os- 
cillators interacting with either two independent reser- 
voirs or a common reservoir was derived perturbatively 
up to the second order of the system-reservoir coupling 
constant (i.e., also in the Born approximation), and the 
entangled dynamics of two-mode Gaussian states was 
analyzed with Ohmic reservoirs as well. In our previ- 
ous work [2^, we have derived the exact master equa- 
tion for two coupled cavity fields under the influence of 
vacuum fluctuation using the Feynman- Vernon influence 
functional theory in the coherent-state path integral for- 
malism poj ). and studied the decoherence dynamics of 
the continuous- variable quantum channel in terms of en- 
tangled two-mode Glauber coherent states with Ohmic 
spectral density. 

In the present work, we shall explore the non- 
Markovian influence of the vacuum fluctuation on the 
continuous-variable quantum channel in terms of an en- 
tangled two-mode squeezed state with different spectral 
densities, i.e. the Ohmic, the sub-Ohmic, and the super- 
Ohmic cases. To study the non-Markovian entanglement 
dynamics of the squeezed-state quantum channel under 
the influence of the vacuum fluctuation, we model the 
system as two cavity fields coupling to a common bosonic 
environment in the at zero temperature. We then use 
the Feynman- Vernon influence functional theory in the 
coherent-state path integral formalism that we have pro- 
vided in our previous work [25| to study nonperturba- 
tively the noise effect on the entanglement dynamics of 
the squeezed states. As is well-known the Feynman- 
Vernon influence functional theory enables us to treat 
both of the back actions from the environment to the sys- 
tem and the system to the environment self-consistently. 
The dissipation and noise dynamics of the quantum chan- 
nel, going beyond the Born-Markov approximation, is 
then governed by an effective action associated with the 
influence functional containing all the influences of the 
environment on the system. 

Utilizing this nonperturbative treatment, the resulting 
exact master equation can be expressed in an operator 
form with time-dependent coefficients describing the full 
dynamics of the back action between the system and the 
environment. We thereby investigate the non-Markovian 
entanglement dynamics of the quantum channel under 
the influence of environments with different spectral den- 



sities, i.e. Ohmic, sub-Ohmic, and super-Ohmic densi- 
ties. The infiuence of the environment induces a shifted 
frequency il{t) and a decay rate r{t) in each cavity mode, 
as well as a shifted coupling strength 17' (t) and a cor- 
related decay rate r'(t) between the two modes. The 
entanglement dynamics depends sensitively on the dif- 
ferent shifted coupling strength Q'(t) (noise- induced en- 
tanglement oscillation), and the decay rates T(t) and 
T'{t) (dissipation-induced suppression of quantum en- 
tanglement) for different spectral densities. We find 
that the Ohmic environment shows a weak dissipation- 
noise effect, while the sub-Ohmic environment leads to 
fast decoherence in the entanglement dynamics. The 
super-Ohmic environment has the strongest memory ef- 
fect, which heavily suppresses the entanglement of the 
squeezed-state quantum channel. 

The paper is organized as follows. In Sec. II, we intro- 
duce the model describing the non-Markovian entangle- 
ment dynamics of the continuous- variable quantum chan- 
nel in terms of an entangled squeezed state, and we shall 
also briefly review the master equation we derived in [25| . 
In Sec. Ill, we use logarithmic negativity as an entan- 
glement measure of continuous- variable states to discuss 
the entanglement dynamics of the squeezed state. The 
numerical results of the entanglement dynamics are given 
in Sec. IV, where we also analyze in detail the influences 
of the environment with different spectral densities on 
the quantum channel. Finally, a brief summary is made 
in Sec. VI. 



II. THE HAMILTONIAN AND THE EXACT 
NON-MARKOVIAN MASTER EQUATION 

A. The model Hamiltonian 

Our system consists of two coupled cavity fields subject 
to a common environment. The Hamiltonian of the total 
system is given by [2^, Hlj 



H = Hs + He + Hi, (1) 

Hs = hujia\ai + fiL02Ci\a,2 + hK{a\a2 + a\ai), 
He = ^ hujkblbk, 

k 

Hi = Y^ Hgikajbt + g*ikaih\), 

Lk 



where 



are the Hamiltonians of the two cavity fields, the environ- 
ment, and their interaction, respectively. The operators 
ai and {I = 1,2) are the corresponding annihilation 
and creation operators of the Ith cavity field with fre- 
quency LOi, and K is a real coupling constant between the 
two cavity fields, which can be realized by a beam split- 
ter. The environment is modeled, as usual, by a set of 
harmonic oscillators described by the annihilation and 
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creation operators bk and = 1, 2, • ■ • ). The couphng 
constants between the cavity fields and the environment 
are given by gik- In the present work we shall consider 
the entangled squeezed state used in the quantum telo- 
protation of continuous- variable states ^14] where the two 
cavity fields are identical, i.e., coi — lu2 ^ luq. We also 
assume that the dominant dissipation and noise effects 
are induced by the vacuum fluctuation so that the en- 
vironment is at zero temperature and the two cavity 
fields should interact homogeneously with the environ- 
ment, namely, gik = g2k = Qk- 

In order to investigate the decoherence effect in the 
entanglement dynamics induced by the environment, a 
specification of the spectral density J(w) of the environ- 
ment is required. The spectral density characterizing the 
coupling strength of the environment to the cavity fields 
with respect to its frequencies is defined by 



tion J{<y.f,<y.'f]t\cx.i,(y.[]Q)^ has the form 



J{^) = ^ \9k?5{<^ - ^k)- 



(2) 



In the continuum limit the spectral density may have the 
form 



JH=,7c^(-)"-iexp(--), 

UJc U}c 



(3) 



where Wc is a cutoff frequency, and "q a dimensionless cou- 
pling constant. The environment is classified as Ohmic 
if n = f, sub-Ohmic if < n < 1, and super-Ohmic if 
n > 1 [4^. 



B. The exact master equation 

The exact master equation describing decoherence dy- 
namics of the two cavity fields can be derived with the 
Feynman- Vernon influence functional method [H, 
in the coherent-state representation [i^. The detailed 
derivation can be found in (25j . and we give only a few 
key steps here for completeness. Going from the quan- 
tum mechanical equation ihdptot (t) /dt ~ [H, ptot (t)] , the 
reduced density matrix fully describing the dynamics of 
the two cavity fields is obtained by integrating out com- 
pletely the environmental degrees of freedom. 



p{a.f,a.'j:;t) — / dfi(cy.i)dfi{cy.'j)j{o!.f,a'j:;t\a.i,a^;0) 



(4) 



where the reduced density matrix is obtained from the to- 
tal density matrix tracing over the environmental degrees 
of freedom, p{af,ay,t) = / dp{z){af,z\ptot{t)\a'j:,z), 
the complex variables a = [ai, 02) and z = (zi, Z2, ■ ■ •) 
the corresponding eigenvalues of the cavity-field opera- 
tors ai_2 and the environment operators bk{k = I, 2, • • ■ ), 
acting on the bosonic coherent-state \a,z), and a de- 
notes the complex conjugate of a. The propagating func- 



J{af,a'j. ;t\ai,a'i;0)^ J D^aD^a.' exp{-{Ss[a, a] 

where S'5[q:, a] is the action of the two cavity fields, 
Ss[a,a] = ^y^i - iaifai{t) + / dT[iai{T)ai{T) 

- A/a/ (r)a/ (r) - Kai {T)av (r)] | , 

and J-[a.,OL,a.' ,a.'] is the Feynman- Vernon influence 
functional obtained after integrating out all the degrees 
of freedom of the environment. 



a, a', a'] = exp I J dr J dr'^ (aj 



a/)(r) 



xp{r - r')«„,(r') + (a/ - a\)ir)p*ir - r')aUr')] }■ 

The time-dependent function /i(r) is the dissipation- noise 
kernel characterizing the full influence of the environment 
on the two cavity fields. 



(6) 



and is completely determined by the spectral density 
J{u})- We only shall consider the spectral density given 
by ([3]) in this paper. 

As we see, all the effects of the environment on the sys- 
tem are incorporated into the influence functional, which 
effectively modifies the action of the cavity system. Since 
the resulting effective action is bilinear in terms of the 
cavity field variables a and ct', the evaluation of the path 
integral over a and a' can be exactly executed with the 
saddle point method. This leads to the dissipation-noise 
equations (I ^ V ) 



ai -I- i(w/a/ -I- nail) 



a[ — i{ijJia[ -\- Ka[i) = — 



2 

r dT'Y,^Ji[r-r' 

2 

rdr'J2^^* ir~r')a'^ (r') , 

•^0 m=l 



(7) 



obeying the boundary conditions ai (0) = an and 
a'l (0) = ajj. The integro-differential dissipation- noise 
equations render the reduced dynamics non-Markovian, 
with the memory of the system interacting with the 
environment registered in the dissipation-noise kernel 
/i(r — r'). Introducing the new variables u{t) and v{t) 
by 

ai (r) = aiiU (r) - ai'iV (r) , 

a[{T) = a[,uiT)-a[,,v{T), I ^ I' , (8) 
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we obtain an explicit solution for the propagating func- 
tion, 

J{af,a'j.;t\ai,a.'^;0) = 

2 

exp I [uaifaii + ua'^a'ij — {uu + vv — l)aJjQ!H] 
1=1 

- ^ [valfav^ + va'na'i, ^ - [uv + vu)a'navi] |. 



(9) 



The non-Markovian master equation can be deduced 
from (gl) and jS]). The resuh is ^ 

+r{t)[2aip{t)a\ — a\aip{t) — p{t)a[ai] 
+r{t)[2a2p{t)al - ala2p{t) - p{t)ala2] 
+r'{t)[2aip{t)al - a\a2pit) - pit)a\a2] 
+r'{t)[2a2pit)a\ - alaip{t) - p{t)alai],{10) 

where 

H'{t) = hn{t){a\ai + 4a2) + hn'{t){a\a2 + alai), 
with 



uu — vv 



n 



- EE -r'(0 - 

(11) 

This is the exact master equation for the dynamics of the 
two cavity fields, fl{t) plays the role of a shifted time- 
dependent frequency for each cavity field, n'{t) accounts 
for a shifted time-dependent coherent coupling between 
the two cavity fields, T(t) represents a time-dependent 
individual decay rate of each cavity field, and T'{t) is a 
correlated decay rate between the two cavity fields. From 
Eq. (fTU)l . we can see that besides the spontaneous decay 
of the individual cavity field, the environment, even if 
only the vacuum fluctuation is concerned, also induces a 
coherent coupling and a correlated spontaneous decay be- 
tween the two cavity fields. The non-Markovian charac- 
ter thus resides in these time-dependent coefficients in the 
master equation. We must emphasize that our derivation 
of the master equation, Eq. ([TO)) , is fully nonperturbative, 
which goes beyond the Born approximation [2^ . [26j and 
involves all the back actions between the environment 
and the cavity fields. 



III. ENTANGLEMENT MEASURE OF 
CONTINUOUS- VARIABLE QUANTUM 
CHANNELS AND ITS DYNAMICS 

A. Logarithmic negativity as entanglement 
measure 

In what follows, we shall analyze the effects of the dif- 
ferent types of noise on the entanglement dynamics of 



the quantum channel in terms of an entangled two-mode 
squeezed state. The entangled two-mode squeezed state 
is defined as the vacuum state acted on by the two-mode 
squeezing operator 



1^(0)) 



_ ^r(aia2- 



-alal) 



|00), 



(12) 



where r is the squeezing parameter. The state approaches 
the ideal Einstein-Podolsky-Rosen (EPR) state in the 
limit of infinite squeezing (r — > oo) [la |. After generat- 
ing the entangled state given by Eq. p2)) . the two cavity 
fields are then propagated, respectively, to the two loca- 
tions separated between the sender and the receiver. The 
quantum channel is thus established through the entan- 
gled two-mode squeezed state and is ready for teleporting 
unknown optical coherent states [ill, [13] ■ The traditional 
way to generate the entangled two-mode squeezed state 
is via the nonlinear optical process of parametric down- 
conversion [l^ . Recently, a microwave cavity QED-based 
scheme to generate such states has also been proposed 

To investigate the entanglement dynamics of the quan- 
tum channel in terms of the two-mode squeezed state, a 
computable entanglement measure for such continuous- 
variable states must be defined first. Here we shall use 
the logarithmic negativity [i^ to quantify the degree of 
entanglement in the quantum channel. The logarithmic 
negativity of a bipartite system was introduced originally 



E 



N 



^Og2^\>^ 



(13) 



where A,^ is the negative eigenvalue of p'^', and p"^' is a 
partial transpose of the bipartite state p with respect to 
the degrees of freedom of the ith party. This measure 
is based on the Peres-Horodecki criterion |46i, |47| that 
a bipartite quantum state is separable if and only if its 
partially transposed state is still positive. 

For the continuous-variable (Gaussian-type) bipartite 
state, its density matrix is characterized by the co- 
variance matrix defined as the second moments of the 
quadrature vector X = {xi,pi,X2,P2), 



(14) 



where AXi = Xi ~ {Xi}, and Xi = 
The canonical commutation relations take the form as 



Pi 



\X^,X^ = iUu, with U ^ 



J \ , . 1 
Q j and J= (^_^ Q 

defining the symplectic structure of the system. The 
property of the covariance matrix V is fully determined 
by its symplectic spectrum ly — [vi, i>2), with ±Vi {vi > 0) 
the eigenvalues of the matrix: iUV. The uncert ainty 
principle exerts a constraint on Vi such that ^ [49|. 
Thus the Peres-Horodecki criterion for the continuous- 
variable state can be rephrased as the state being sepa- 
rable if and only if the uncertainty principle, V+^U ^ 0, 
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is still obeyed by the covariance matrix under the partial 
transposition with respect to the degrees of freedom of 
a specific subsystem [4a|- In terms of phase space, the 
action of partial transposition amounts to a mirror re- 
flection with respect to one of the canonical variables of 
the related subsystem. For instance, V = AVA , and 
A = diag{l, 1, 1, —1) is the partial transposition with re- 
spect to the second subsystem. If a Gaussian-type bipar- 
tite state is nonseparable, the covariance matrix V will 
violate the uncertainty principle and its symplectic spec- 
trum D — 1/2) will fail to satisfy the constraint ^■ 
The logarithmic negativity is then used to quantify this 
violation as [i^ 

En = max{0, - log2(2j>„iin)}, (15) 

where v-adn is the smaller one of the two symplectic eigen- 
values Vi. It is evident from Eq. that, if V obeys 
the uncertainty principle, i.e., Vi ^ i, then E]y{p) = 0, 
namely, the state is separable. Otherwise, it is entan- 
gled. Therefore, the symplectic eigenvalue Vmin encodes 
a qualitative feature of the entanglement for an arbitrary 
continuous-variable bipartite state. 



B. The entanglement dynamics 

With this entanglement measure at hand, we study 
now the entanglement dynamics of the squeezed-state 
quantum channel in our model. A straightforward way 
to obtain the time-dependent solution of the entangled 
squeezed state is by integrating the propagator function 
over the initial state of Eq. (|4]) , where the initial state in 
coherent-state representation is given by 

, exp[- tanhr(ai,a2i + a'iia'2i)\ 

p{OL^, a- 0) = -2 . (16) 

cosh r 

The solution of the reduced density matrix can be ob- 
tained exactly, 

p[a.f,a.'f;t) = 69 cxp[^(6ia^j + 62";/ + hai/a'if 
i^i' 

+ b5aifa'i,f + —aifauf + — aj^aj,^)], (17) 

where the time-dependent parameters bi {i — 0, ■ ■ ■ ,6) 
are given explicitly in the Appendix. 

From the above solution, the covariance matrix V 
can be calculated analytically, and the logarithmic neg- 
ativity E^it) can also be obtained exactly from Eq. 
(fT5|) . It is easy to verify that the initial entanglement 
is En (0) — . While the asymptotical behavior of the 
entanglement (in the long time limit) can be found from 
the solution 

p{t ^ Co) = \'lpasy.){'ipasy.\, 

li'asy.) = V&[^e'''i('^5-'^^)>0), (18) 



where b'^ = bf,{t -> 00) = -^^^ and b\ = bi{t ^ 00) = 
tanhre ^ ' » — ^_ ^ This asymptotical solution results in 

i?Ar(i^(^) = ^, (19) 

namely, the final entanglement is only one half of the ini- 
tial value. Equation indicates that the environment 
does decrease the entanglement of the quantum channel, 
as a dissipation-noise effect. This asymptotical result is 
also consistent with the Markovian limit at zero temper- 
ature in since the non-Markovian dynamics must be 
asymptotically reduced to the Markovian limit |25| . 

It should be noted that the noise behavior of the quan- 
tum channel also depends on the structure of the ini- 
tial two-mode squeezed state. This can be easily seen if 
we introduce the operators related to the center-of-mass 
and relative motional variables of the two cavity fields as 
= {a\ + al) and = {a\ — al), respectively. Then 
the initial state can be rewritten in terms of these two 
operators, 

IV'(0)) = ^-e-^(^^^-'^^^)|00). (20) 
coshr 

Since the two cavity fields interact homogeneously with 
the environment, namely, gik — g2k — gk, the interaction 
between the cavity fields and the environment only influ- 
ences the dynamics of the center-of-mass variable; it has 
no effect on the relative motion of the two cavity flelds (a 
similar discussion for two harmonic oscillators interacting 
with a common reservoir is given in [131 )■ other words, 
the part of the squeezed state relating to the relative 
variable is immune to the environment, while that of the 
center-of-mass experiences severe dissipation and noise 
from the environment. This results in the solution (jlSp . 
One can verify that Eq. is indeed a decoherence-free 
squeezed state, i.e., \^dfs) c< e^'^'^i"'^^) |00). If such a 
state serves as the quantum channel, the quantum chan- 
nel is free from the vacuum fluctuation, and the entan- 
glement is preserved. 

IV. NUMERICAL ANALYSIS OF THE 
NON-MARKOVIAN ENTANGLEMENT 
DYNAMICS 

The full entanglement dynamics of the squeezed-state 
quantum channel is determined by the reduced density 
matrix obeying the master equation (jlOp . To solve the 
master equation, we must flnd first the time-dependent 
coefficients contained in the master equation, the shifted 
frequency ri(i) and the shifted coherent coupling fl'{t), 
as well as the individual and correlated decay rates T{t) 
and r'(t). These coefficients are completely determined 
by the functions u{t) and v{t) as the solutions of the 
dissipation- noise equations ([7]) via ([5]). However, the 
dissipation-noise equations have to be solved numerically 
for the general environmental spectral density ([3]). 



6 



0.4 r 

I 

0.3 ' 
§ 0.2 i 




0.2 0.4 0.6 0.8 1 

FIG. 1: Time-dependence of the frequency shift S^l{t) {— 
SQ'{t)) induced by the environment with spectral density (|3]) 
of n = 1 for the Ohmic (solid line), n — 3 for the super- 
Ohmic (dashed line), and n — 1/2 for the sub-Ohmic (dotted 
line). The asymptotical value 5Q,{t oo) — rjujc, \/ttt}ujc and 
2riu)c for the Ohmic, the super-Ohmic, and the sub-Ohmic, 
respectively. The parameters in Q are taken as r; = 0.005 
and ujc = SO.Oli^o, while the coupling constant between the 
cavity fields as k = O.Sojq. 



In Figs. [T]and[H wc plot the numerical results for the 
frequency shift Sft{t) = luq — ^l{t) and decay rate T{t) 
of the individual cavity field. We choose three different 
spectral densities: n = 1, 1/2, and 3 for the Ohmic, the 
sub-Ohmic, and the super-Ohmic spectral densities, re- 
spectively. Since the two cavity fields are considered to 
be identical {uji = W2 =1-^0) a-nd interact homogeneously 
with the common environment {gik = g2k = gk), the 
environment-induced shifts of the field frequencies and 
the coherent coupling between the two cavity fields are 
equal, i.e., Sn{t) = 6n'{t) where 6n'{t) = k~ n'{t). The 
individual and correlated decay rates are also equal to 
each other, r{t) = T'{t), for the same reason. From Figs. 
[T] and [51 we find that the dissipation-noise dynamics is 
characterized by two time scales: ri = (the short- 
est time scale of the environment) and T2 = oJq^ (the 
time scale of the cavity fields). When t < ti, both coef- 
ficients (5r2(i) and T{t) grow very quickly. After ti, they 
approach the corresponding asymptotical values gradu- 
ally in the time scale of T2. We should point out that 
the asymptotic values of Sn{t) and T{t) in the Ohmic 
spectral density reproduce the Markovian limit, as we 
have shown in our previous work psj). Compared with 
the super-Ohmic case, the sub-Ohmic dissipation shows 
a slower asymptotical tendency in the time scale T2 . This 
is because, in the super-Ohmic case, the short-time cor- 
relation (the ultraviolet mode) is dominant, while in the 
sub-Ohmic case, the long-time correlation (the infrared 
mode) becomes important. 

In fact. Fig. [2] also tells us that the decay rate 
grows very fast in the time scale ti and develops a jolt 
[39j . This peak manifests the significant effect of the 
non-Markovian dynamics, especially for the super-Ohmic 
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FIG. 2: Time-dependence of the decay rate r{t) (= r'{t)) 
with n — 1 for the Ohmic (solid line), n — 3 for the super- 
Ohmic (dashed line), and n = 1/2 for the sub-Ohmic (dotted 
line) environments. The input parameters in the numerical 
calculation are the same as in Fig. 1. 

case. In the usual Born-Markovian approximation used 
in the literature [13, [H, [13, [IH , the back action of the 
environment on the system is completely ignored by the 
assumption of the response time of the environment be- 
ing much smaller than the characteristic time (T2) of the 
system. The decay rate becomes then time independent. 
The time dependence of the decay rate in the exact mas- 
ter equation PH)) contains the full memory effect of the 
system interacting with the environment, as a result of 
the non-Markovian dynamics. The non-Markovian en- 
tanglement dynamics of the quantum channel thus be- 
comes transparent due to the presence of the time de- 
pendence of these coefficients in the time scale ti. This 
short-time correlation will infiuence strongly the later- 
time entanglement dynamics of the squeezed state. 

In Figs. [3] and m we plot the time evolution of the loga- 
rithmic negativity for the entanglement dynamics of the 
squeezed state. Figure [3] describes the case where the 
two cavity fields are initially decoupled (k = 0). As one 
would expect, in the absence of vacuum fluctuations the 
entanglement does not change in time. Figure 2] shows 
the coupling cavity fields (k = 0.5), where the entangle- 
ment undergoes a lossless periodic oscillation. When the 
vacuum fiuctuation is taken into account, the entangle- 
ment dynamics is significantly changed as we see from 
Figs. [3] and m On one hand, the vacuum fluctuation in- 
duces entanglement oscillations (see Fig. ^ or shifts the 
frequencies of the entanglement oscillations (see Fig. 
due to the effect of the shifted two- mode coupling Sil'{t) 
(given in Fig. [T]). We should mention that the entan- 
glement oscillations have also been found for two uncou- 
pled harmonic oscillators interacting with two indepen- 
dent reservoirs On the other hand, the amplitude 
of the entanglement oscillation is suppressed gradually 
and tends to its asymptotical value (= r/ln2) arising 
from the joint dissipation effects of the individual and 
correlated decay rates T{t) and r'(t). Furthermore, the 
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5 10 15 20 25 

FIG. 3: Time evolution of the logarithmic negativity EN{t) 
without the noise (dash-dotted line), and with the Ohmic 
(solid line), the super-Ohmic (dashed line), and the sub- 
Ohmic (dotted line) noise environments, in which the two 
cavity fields initially have no coupling to each other (k = 0). 
The other input parameters are still the same, and the squeez- 
ing parameter r — 3.0 
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FIG. 4: Time evolution of the logarithmic negativity i?jv(f) 
without the noise effect (dash-dotted line), and with the 
Ohmic (solid line), super-Ohmic (dashed line), and sub- 
Ohmic (dotted line) noise environments. The input parame- 
ters are the same as in Figs. 3, except for k — 0.5a;o 



dissipation will also erase the oscillation of entanglement 
arising from the coherent coupling between the two cav- 
ity fields. These decoherence effects are consistent with 
that obtained from the discrete qubit models [5l|. It 
shows again that the asymptotic value of the entangle- 
ment reproduces the result in the Markovian limit ^17|. 

One may also see from Figs. [3] and 2] that the en- 
tanglement dynamics behaves differently for three differ- 
ent spectral densities, as well as for the two coupled and 
uncoupled cavity fields. Comparing with the oscillation 
behaviors presented in Figs. [3] and 01 it shows that the 
order of the entanglement oscillations with three differ- 
ent spectral densities is reversed for the coupled and un- 



coupled cavity fields. This behavior comes mainly from 
the fact that the super-Ohmic environment induces the 
strongest entanglement oscillation; next is the sub-Ohmic 
case; while the Ohmic environment causes a relatively 
weak entanglement oscillation, as shown in Fig. [T] By 
the definition n'{t) = k — Sn'{t), it is easy to check that, 
for the uncoupled cavity fields, k = so that |f^stip(^)l > 
I^sm6(0I > \^'ohmi^)\- This leads to the frequencies of the 
entanglement oscillations fsup > fsub > fohm, as shown 
in Fig. [21 While for the coupled cavity fields, k = 0.5 
so that \n',^^{t)\ < < |f^;,„(t)|. This gives the 

frequency ordering fg^p < fsub < fohm, as plotted in 
Fig. m This indicates that the non-Markovian induced 
entanglement oscillations can be quite different for the 
two coupled or uncoupled two entangled cavity fields. 

Meanwhile, the entanglement oscillation is sustained 
for the longest time in the Ohmic case, while the super- 
and sub-Ohmic environments cause severe decoherence. 
The sub-Ohmic (low-frequency) vacuum fluctuation in- 
duces a strong dissipative dynamics (the largest decay 
rate, as in Fig. [J), and therefore results in a fast de- 
coherenc of the entanglement of the squeezed state, as 
expected. However, a remarkable result occurs in the 
super-Ohmic case where the entanglement and its oscil- 
lation are also strongly suppressed. In contrast to the 
sub-Ohmic case, the decay rate in the super-Ohmic case 
is almost negligible except for a sharp peak in the short- 
time scale Ti, as shown in Fig. [2l We find that this 
short-time sharp peak induces a significant contribution 
to the entanglement decoherence in the quantum chan- 
nel. This decoherence effect is a manifestation of the 
memory dynamics between the system and the environ- 
ment. It is this non-Markovian dynamics that causes a 
rapid decoherence of the entanglement of the squeezed 
state in the super-Ohmic environment. We should note 
that, within the time scale ti = lo^^, the initial peak of 
the decay rates also exists in the Ohmic and sub-Ohmic 
cases although it is not as strong as in the super-Ohmic 
case. This initial "jolt" in the decay rates is indeed a 
general feature of the non-Markovian processes for the 
decoherence enhancement, as pointed out first by Hu et 
al. [39] in the study of quantum Brownian motion. The 
results we obtained in this work demonstrate that the 
non-Markovian dynamics also speed the decoherence of 
the entanglement in continuous-variable quantum chan- 
nels. 



V. SUMMARY 

In the present work, we have studied the detrimen- 
tal effects of the environment on the continuous- variable 
quantum channel in terms of the entangled two-mode 
squeezed state. Using the Feynman- Vernon influence 
functional theory in the coherent-state path integral rep- 
resentation, we derive the exact master equation for the 
two cavity fields under the influence of vacuum fluctua- 
tion psj and then investigate the non-Markovian entan- 
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glement dynamics of the two-mode squeezed state quan- 
tum channel utihzed in quantum teleportation [l4{ for 
three different spectral densities, the Ohmic, the sub- 
Ohmic, and the super-Ohmic non-Markovian environ- 
ments. Very recently, a similar exact master equation 
has also been derived for two harmonic oscillators lin- 
early coupled to a thermal bath where the entanglement 
dynamics is studied in the Markovian approximation [2^ . 

We numerically study the non-Markovian entangle- 
ment dynamics of the quantum channel based on the 
exact master equation ([TUl) for three different noise en- 
vironments. Our numerical result indicates that the en- 
tanglement dynamics behaves different for the different 
environmental spectral densities which leads to signifi- 
cant distinctness in the time-dependent behavior of the 
dissipation-noise function, in particular, within the short 
time scale ri of the environment. For Ohmic environment 
the system shows the longest quantum coherence because 
of the weak time-dependent dissipation of the entangle- 
ment dynamics of the squeezed state. In the sub-Ohmic 
case the squeezed state has a strong dissipation dynamics 
(corresponding to a large decay rate) induced mainly by 
the low- frequency noise of the environment, which results 
in fast decoherence for the entanglement dynamics of the 
squeezed state. The most significant evidence of the non- 
Markovian dynamics occurs in the super-Ohmic environ- 
ment in which the strong non-Markovian process near the 
short time scale (ri) speeds the decoherence of the en- 
tanglement. These non-Markovian properties are indeed 
consistent with the non-Markovian phenomena explored 
in quantum Brownian motion [39| . 

We may also point out that for the squeezed-state 
quantum channel considered in this paper, both the 
asymptotical and the numerical solutions show that one 
half of the initial entanglement carried by the squeezed 
state will be retained regardless of the spectral density of 
the environment. This is consistent with the solution in 
the Markovian limit ^17.] . This result depends only on the 



structure of the initial squeezed state as well as the prop- 
erty of the homogeneous coupling between the system 
and the environment. Thus the asymptotical state 
is indeed a decoherence-free entangled squeezed state 
in our model, which may serve as a noiseless quantum 
channel for further applications in quantum communi- 
cation. But it should be pointed out that, if the two 
cavity fields couple with two independent reservoirs, the 
above decoherence-free state no longer exists and the re- 
maining entanglement will eventually be lost completely 
[1^. As our concentration is on the optical cavity fields, 
we have considered only the zero-temperature environ- 
ment. A more general case, e.g., with the environment 
at a finite temperature, could hopefully be figured out 
by a similar approach to the derivation of Eq. ©. As 
robustness of the quantum channel is essential in view 
of decoherence, we hope that our consideration of non- 
Markovian entanglement dynamics in this paper provides 
useful information for experimental designs of quantum- 
information protocols. 
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Appendix: the coefficients of p{a.f,ay,t) and their 
asymptotical behaviors 

The explicit form of p{(y f , a'y ; t) is obtained from Eq. 
(ID) by the evaluation of the integration. The final solu- 
tion is given by Eq. (|f 7p with the parameters 
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cosh^ 1 — 2 tanh^ r(m^ + n^) + tanh^ r(rn^ — n^)^ 



etanh r (un + wm) + (um + wn) + c tanh r{un + vm){um + vn) 

Oi = 2 4 ^ -+iBxi!iruv, 

1 — 2 tanh r{m? + n^) + tanh r(m^ — n^)^ 
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etanh^ r{v? + w^) + c tanh rut; 



1 — 2 tanh^ r{rnP' + n^) + tanh"* r{m? —n^y 
r(un + vm){um + vn)] + c{— tan 
1-2 tanh^ r{m?' + rfi) + tai 
e{— 2tanh^ r[u{un + vm) + v{um + vn)]} + c{— tanh^ r[u{um + vn) + v{un + um)]} 

1 — 2 tanh^ r{m? + n^) + tanh^ r{vn? — n^)^ 
e{2tanh'' r[u{um + vn) + v{un + vm)]} + c{tanh^ r[u{un + vm) + v{um + vn)]} 
1 — 2 tanh^ r(m^ + n^) + tanh^ r(m^ — n^)^ ' 
e(— 4tanh^ ruu) + c[— tanhr(u^ + 1)^)] 



e[— 4tanh r{un + vm){um + vn)] + c{— tanh r[{un + vm) + {um + vn) ]} , / 2 2\ 

03 = 2 — 7 — o o\ 4 — 7 — o oTo tannr(M +v ), 

1 — 2 tanh r(m^ + n^) + tanh r{m^ — n^)^ 

o S O "l- Q n Vi ^ IT' l/Ti ( • 
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1 - 2 tanh^ r(m2 + n^) + tanh* r(m2 - n^)^ ' 



where c = 1 — tanh r(m^ + n^), e = tanh rmn, m = 
UU + vv — 1, and n = uv -\- vu. In the long time Hmit, 
u{t — »• 00) = v{t — »• 00) = ^ Then 

6o(i ^ 00) = — \ — , bAt — > 00) = b^it — > 00) = 0, 
coshr 

Oi(f— >oo) = 02(4^00) = — 

-b^{t 00) tanhre-2»('^')-'=)* 
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